Linear Response, Validity of Semi-Classical Gravity, 
and the Stability of Flat Space 



Paul R. Anderson^ *, Carmen Molina-Paris^''^ and Emil Mottola^ ^ 
^ Department of Physics, Wake Forest University, Winston-Salem, North Carolina, 27109 
T-8, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545 
Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK 
(LA-UR-02-5809) 



(N 



A quantitative test for the validity of the semi-classical approximation in gravity is given. The 
criterion proposed is that solutions to the semi- classical Einstein equations should be stable to lin- 
earized perturbations, in the sense that no gauge invariant perturbation should become unbounded 
in time. A self-consistent linear response analysis of these perturbations, based upon an invariant 
effective action principle, necessarily involves metric fluctuations about the mean semi-classical ge- 
^— ^ ■ ometry, and brings in the two-point correlation function of the quantum energy-momentum tensor 

' in a natural way. This linear response equation contains no state dependent divergences and requires 

, no new renormalization counterterms beyond those required in the leading order semi-classical ap- 

_ ' proximation. The general linear response criterion is applied to the specific example of a scalar 

field with arbitrary mass and curvature coupling in the vacuum state of Minkowski spacetime. The 
, spectral representation of the vacuum polarization function is computed in n dimensional Minkowski 

spacetime, and used to show that the flat space solution to the semi-classical Einstein equations for 
n = 4 is stable to all perturbations on distance scales much larger than the Planck length. 
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I. INTRODUCTION 



> 

in 

There are many well known difficulties that arise when attempting to combine quantum field 
theory and general relativity into a full quantum theory of gravity. Almost certainly, a consistent 
quantum theory at the Planck scale requires a fundamentally different set of principles from those 
of classical general relativity, in which even the concept of spacetime itself is likely to be radically 
altered. Yet, over a very wide range of distance scales, from that of the electroweak interactions 
(10~^^ cm) to cosmology (10^^ cm), the basic framework of a spacetime metric theory obeying 
^ ■ general coordinate invariance is assumed to be valid, and receives phenomenological support both 

from the successes of flat space quantum fleld theory at the lower end of this distance scale, and 
5-H ■ classical general relativity at its upper end. Hence, whatever the full quantum theory of gravity 

entails, it should reduce to an effective low energy fleld theory on this very broad range of some 43 
orders of magnitude of distance [1,50]. 
. , ■ To the extent that quantum effects are relevant at all in gravitational phenomena within this 

' range of scales, one would expect to be able to apply semi-classical techniques to the low energy 

effective theory of gravity. In the semi-classical approximation to gravity the spacetime metric 
gab is treated as a classical c-number fleld and its quantum fluctuations are neglected, although 
quantum fluctuations of the other flelds are taken into account. The semi-classical approach has 
been discussed and studied for some time now, and a considerable body of results has been obtained 
[3]. Yet a definitive answer to the question of what is the limit of validity of this approach, and 
under what physical conditions it must break down has remained somewhat unclear. 

It is our purpose in this paper to propose a well-defined quantitative criterion for the validity of 
the semi-classical approximation to gravity, within the semi-classical formalism itself, namely that 
solutions to the semi-classical Einstein equations should be stable against linearized perturbations 
of the geometry. This criterion may be formulated by means of a linear response analysis [4-6]. 

It is important to distinguish what we mean in this paper by the semi-classical approximation 
from the ordinary loop expansion, which is sometimes also called semi-classical. In the ordinary 
loop expansion of the effective action, h is the formal (loop expansion) parameter. As a result both 
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the matter and gravitational quantum fluctuations are treated on exactly the same footing, and the 
back-reaction of these fluctuations on the metric (being first order in the expansion parameter h) are 
neglected. If one does attempt to include such effects in some modified loop expansion, the technical 
issues involved in defining a one-loop effective action for gravitons that respects both linearized 
gauge and background field coordinate invariance must be faced. These are difficult enough to have 
impeded progress in the standard loop expansion in gravity [7,8]. An unambiguous definition of the 
corresponding conserved and gauge invariant energy-momentum tensor for gravitons in an arbitrary 
curved spacetime has not yet been given [9]. Apart from such technical difficulties it should be 
clear that a simple loop expansion is ill-suited physically to many applications that have been and 
are likely to be of interest in semi-classical gravity, such as particle creation in the early universe, 
or black hole radiance, where the quantum effects of matter significantly affect the background 
geometry after some period of time, but (it is usually assumed) the quantum fluctuations of the 
geometry itself can be neglected. Whenever quantum effects of matter are expected significantly to 
affect the classical geometry, the standard loop expansion, which treats these effects as order h and 
small, must certainly break down. 

The semi-classical approximation to gravity wc discuss in this paper treats the matter fields as 
quantum but the spacetime metric as classical, and allows for the consistent back-reaction of the 
quantum matter on the classical geometry. This asymmetric treatment can be justified formally 
by replicating the number of matter fields A'^ times and taking the large N limit of the quantum 
effective action for the matter fields in an arbitrary background metric Qat [10]. Then, the semi- 
classical equations for the metric are derived by varying the effective action, with local gravitational 
terms included. Since no assumption of the weakness or perturbative nature of the metric is assumed, 
the large N expansion is able to address problems in which gravitational effects on the matter are 
strong, and the matter fields can have a large cumulative effect on the classical geometry in turn. The 
absence of quantum gravitational effects in the lowest order large N approximation also means that 
the technical obstacles arising from the quantum fiuctuations of the geometry are avoided. General 
coordinate invariance is assured, provided only that the matter effective action is regularized and 
renormalized in a manner which respects that invariance [7]. In that case the quantum expectation 
value of the matter energy-momentum tensor (Tab) is necessarily conserved. 

Assuming that the classical energy-momentum tensor for the matter field(s) vanishes (an as- 
sumption that may be easily relaxed if necessary) , the unrenormalized semi-classical back- reaction 
equations take the form [11] 

Gab + J^Qab = 87rG„ {Tab) ■ 

Here Gab is the Einstein tensor, A is the cosmological constant (which may be taken to be zero 
in some applications), Gj^, is Newton's constant, and {Tab) is the expectation value of the energy- 
momentum tensor operator of the quantized matter ficld(s) . Among the technical issues that must be 
confronted is the renormalization of the expectation value of Tat, a quartically divergent composite 
operator in n = 4 spacetime dimensions. The renormalization of its expectation value requires the 
introduction of fourth order counterterms in the effective action, that modify the geometric terms 
on the left hand side of eq. (1.1) [3]. 

Once a renormalized semi-classical theory has been defined, one possible route to investigating 
its validity is to compare calculations in a theory of quantum gravity with similar semi-classical 
calculations. Since a well-defined, full quantum theory is lacking, this has been done only in some 
simplified models of quantum gravity. Large quantum gravity effects were found in three-dimensional 
models by Ashtckar [12] and Beetle [13]. In four dimensions. Ford has considered the case of graviton 
production in a linearized theory of quantum gravity on a flat space background, and compared the 
results with the production of gravitational waves in semi-classical gravity [14]. He found that they 
were comparable when the renormalized energy-momentum (connected) correlation function, 

{Tab{x)T,d{x'))co^ = {Tab{x)T,d{x')) - {Tab{x)) {Tcd{x')) 

satisfied the condition 

(T„6(a;)Ted(a;'))co„ <C {Tab{x)){T,d{x')) . 

The limits of validity of the semi-classical approximation have also been studied without making 

reference to a speciflc model of quantum gravity. Kuo and Ford [15] proposed that a measure of 
how strongly the semi-classical approximation is violated can be given by how large the quantity, 

{Tab{x)T,d{x'))c^r, 

{Tab(x)T,d(x')) 



Aabcd{x,x') = 



is, where it is assumed that the expectation values in this expression arc suitably rcnormalized. It is 
important to note that eq. (1.4) is coordinate dependent, since both the numerator and denominator 
are tensor quantities. The situation is complicated further by the regulaxization and renormalization 
issues that arise in defining the quantities appearing in this expression. Using normal ordering, Kuo 
and Ford [15] computed the quantity 

con 

(roo(a;)Too(a;)) 

for a free scalar field in flat space for several states including the Casimir vacuum. They found that 
it vanishes in a coherent state, whereas in many other cases, including the Casimir vacuum, it is of 
order unity. 

Wu and Ford [16] computed the radial flux component of cq. (1.4), in the cases of a moving mirror 
in 2-dimensions and an evaporating black hole far from the event horizon in both 2 and 4 dimensions. 
They found that it was of order unity over time scales comparable to the black hole mass, but that 
it averages to zero over much larger times. In a normal ordering prescription they found state 
dependent divergent terms. They also showed that in the simple case of radiation exerting a force 
on a mirror, the quantum fluctuations in the radiation pressure are due to a state dependent cross 
term in the energy-momentum tensor correlation function [17,18]. 

Phillips and Hu [19] used zeta function regularization to compute A(a;) with the denominator 
replaced by the quantity (roo(a;))^, for a free scalar field in some curved spacetimes having Euclidean 
sections. They also computed A{x) for a scalar field in flat space in the Minkowski vacuum state, 
using both point splitting and a smearing operator to remove the divergences [20]. For the flat 
space calculation they found that A(x) depends on the direction the points axe split, but that it is 
of order unity regardless of how the points are split. They used their results to criticize the Kuo- Ford 
conjecture and to suggest that the criteria for the validity of the semi-classical approximation should 
depend on the scale at which the system is being probed. 

Wu and Ford [18] addressed the Kuo-Ford conjecture and the above mentioned criticism of it 
by Phillips and Hu. They stated that the conjecture is incomplete because it does not address the 
effect of divergent state dependent terms. They suggested that any criterion for the validity of the 
semi-classical approximation should be a non-local one that involves integrals over the world lines 
of test particles. They also argued that the question of whether the semi-classical approximation is 
valid depends on the specifics of a given situation, including the scales being probed and the choice 
of initial quantum state. 

Although it is somewhat unclear what the dimensionless small parameter is that controls the 
inequality (1.3), Ford's initial work and these subsequent discussions draw attention to the impor- 
tance of the higher point correlation functions of the energy-momentum tensor. It is quite clear, at 
least in qualitative terms, that if the higher point connected correlation functions of Tab are large 
(in an appropriate sense to be determined), it cannot be correct to neglect them completely, as the 
semi-classical equations (1.1) certainly do. 

Another context in which the quantity {Tab{'-c)Tcd{x')) con plays a role is stochastic semi-classical 
gravity [1,21-25]. In this case the probability distribution function for the quantum noise is obtained 
from the symmetric part of this correlation function [21]. The dissipation kernel has been shown 
to be related to the antisymmetric part of this correlation function [26]. Stochastic semi-classical 
gravity is an interesting attempt to go beyond the semi-classical approximation. However, for the 
purposes of the present work, we do not make any stochastic assumptions and investigate the validity 
of the semi-classical equations within the large A'' approximation itself. 

The energy-momentum correlation function {Tab{x)Tcd{x'))con has been directly computed for a 
scalar field in a two dimensional spacetime with a moving boundary [27], for scalar fields and the 
Maxwell field in Minkowski spacetime [28,29], and for a masslcss minimally coupled scalar field in de 
Sitter spacetime, in the case that the points are spacelike separated and geodesically connected [30]. 
It has also been computed indirectly through the non-local kernel appearing in the deviation of (Tab) 
from flat space [31-33], from a Robertson- Walker spacetime [34], and from a general conformally 
flat spacetime [35]. The noise and dissipation kernels in stochastic semi-classical gravity are related 
to the energy-momentum tensor correlation function [26]. These quantities have been computed 
exactly or approximately for scalar fields of various types in several situations including Minkowski 
spacetime [36,37], hot flat space [26], the far field limit of a black hole in equilibrium with a thermal 
field [38], Robertson- Walker spacetimes [21,22,39,40], Bianchi type I spacetimes [23], and a weakly 
curved spacetime using a covariant expansion in powers of the curvature [41]. 

Although technical problems such as renormalization and coordinate invariance complicate mat- 
ters, this body of previous work suggests that the correlation function (Tai,(a;)rcd(a;'))con should 



A{x)^ 



play an important role in determining the validity of the semi-classical approximation. However, 
the proper context for incorporating and making use of the information contained in this correlation 
function in a well-defined {i.e., finite), quantitative framework, that respects general coordinate 
invariance has remained somewhat unclear. 

The criterion we propose in this paper, that solutions to the semi-classical Einstein equations 
should be stable against linearized perturbations of the geometry, provides just such a framework. 
According to standard linear response theory [4,5,42], the linearized equations for the perturbed 
metric depend on the retarded two-point correlation function of the energy-momentum tensor eval- 
uated in the semi-classical background metric gab [6]. In this case, the correlation function can be 
computed using the Closed Time Path (CTP) effective action [43]. The result is a retarded cor- 
relation function that involves the commutator of two energy-momentum tensor operators. Hence 
the perturbations are manifestly causal. Moreover, the UV divergences found in the unrenormalized 
linear response equations are exactly those required to renormalize the semi-classical theory itself. 
This ensures that no state-dependent divergences occur. Finally, gauge transformations of the lin- 
earized metric fluctuations, hat, are easily handled within the linear response framework, so that 
ambiguities related to quantities such as (1.4) do not arise. Thus, standard linear response theory 
provides a well-defined framework in which to test the validity of the semi-classical approximation 
to gravity, which directly involves {Tabix)Tcd{x')} and its renormalization, in a manner that is in 
complete accordance with the physical principles of general covariance and causality. 

Since this criterion for the validity of the semi-classical approximation lies strictly within the 
context of that approximation itself, one avoids problems such as gauge invariance of the energy- 
momentum tensor for gravitons, that inevitably appear if one tries to go beyond the semi-classical 
approximation and include quantum effects due to the gravitational field. Although these effects 
certainly are not contained in the semi-classical Einstein eqs. (1.1), it is possible to study the 
properties of linearized gravitational fluctuations about the self-consistent solution of cq. (1.1), 
simply by taking one higher variation of the effective action that leads to that etiuation. This 
second variation involves the two-point correlation function (1.2), evaluated in the self-consistent 
background geometry. 

To understand qualitatively the role of the two-point correlation function in the validity of the 
semi-classical approximation, it is helpful to consider the physical analogy between semi-classical 
gravity and semi-classical electromagnetism. The connected correlation function (1.2) measures 
the gravitational vacuum polarization, which contributes to the proper self-energy of the linearized 
graviton fluctuations around the background metric, just as the current two-point correlation func- 
tion, {j°{x)j''{x'))con, measures the electromagnetic vacuum polarization which contributes to the 
proper self-energy of the photon [44,45]. Hence, if these polarization effects are significant, the 
semi-classical approximation has certainly broken down, at least in the form specified by eq. (1.1), 
where all fluctuations of the metric have been ignored. In quantum electromagnetism (QED) we 
know exactly how to take these fluctuation effects into account, namely by scattering and interac- 
tion Feynmaii diagrams involving the photon propagator. These processes arc important not only in 
scattering between a few particles at high energies, but also in low energy processes in hot or dense 
plasmas [5]. Analogous statements should be applicable to gravity. Thus, if the linear response 
validity criterion is not satisfied, there will be no avoiding the technical difficulties and physical con- 
sequences of treating the fluctuations of the gravitational field itself, even if we seek to understand 
only the infrared behavior of a semi-classical approximation to the effective theory of gravity, far 
below the Planck energy scale. 

As a particular illustration of the validity criterion, we apply it to the example of a scalar field 
with arbitrary mass and curvature coupling in the vacuum state of Minkowski spacetime. We 
express the retarded correlation function of the linear response analysis in fiat space in terms of 
a Kallen-Lehmaun spectral representation [45]. The positivity of the spectral representation is 
sufficient to demonstrate that there are no unstable modes of the linearized semi-classical equations 
around flat space at distance scales far larger than the Planck scale, and hence, that flat spacetime 
is completely infrared stable in semi-classical gravity. The semi-classical stability of Minkowski 
spacetime has been investigated previously by several authors [32,33,46], and instabilities have been 
found which involve strictly Planck scale variations of the metric fluctuations in space and/or time, 
which arise from the terms fourth order in derivatives of the metric that are needed to renormalize 
eqs. (1.1). Their existence clearly precludes the validity of the semi-classical large N approximation 
at Planck time or distance scales. Prescriptions for explicitly reducing the order of the equations, 
which eliminates these Planck scale solutions have been proposed [47,48] and discussed in some 
detail [49]. Whether or not these prescriptions are accepted in the general case, it is quite clear 
a priori that the semi-classical approximation (1.1) can be viewed at best only as the low energy 
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effective field theory limit of a more complete quantum theory [1,50], and that no reliable results 
can be obtained from this approximation in the Planckian regime. However, the flat space example 
treated in some detail in this paper shows explicitly that the semi-classical approximation does give 
mathematically meaningful and physically sensible results, when properly restricted to its range of 
validity at space and distance scales very much larger than the Planck scale. 

The organization of the paper is as follows. In the next Section the properties of the large N 
semi-classical approximation in gravity and its renormalization within the covariant effective action 
framework are reviewed. In Section III the linear response theory for the semi-classical back-reaction 
equations is described. The form of the two-point correlation function for the energy-momentum 
tensor that appears in the linear response equations is given, and its properties and renormalization 
arc discussed. Then our proposal for a necessary condition for the validity of the semi-classical 
approximation is presented. In Section IV the use of our criterion is illustrated for the case of a 
scalar field with arbitrary mass and curvature coupling in the vacuum state in Minkowski spacetime. 
The linear response analysis implies that flat spacetime is stable under small fluctuations at large 
wavelengths. Our results are discussed further in Section V. Some additional applications of our 
criterion to the study of quantum effects in cosmological and black hole spacetimcs arc suggested. 
There are two Appendices. The first deals with the general decomposition of tensors and polarization 
operators in Minkowski spacetime. The second contains the technical details of the computation of 
the retarded correlation response function for a scalar field in Minkowski spacetime. 

II. SEMI-CLASSICAL GRAVITY AND RENORMALIZATION 

The most direct route to the semi- classical equations (1.1) is via the effective action method in the 
large N limit. We consider the specific example of N non-interacting scalar fields. Generalizations 
to interacting fields and fields of other spin are straightforward, but as they are not required to 
expose the main elements of the stability criterion, we treat only this simplest case in detail. We 
begin by reviewing the effective action formulation of the semi-classical eqs. (1.1) without regard to 
boundary conditions or the state of the field. Thus, the equations in this section are valid for both 
the (out\in) and (m|m) formalisms. We postpone to the next Section the introduction of the OTP 
method which selects real and causal {in\in} expectation values. It is this latter form that must be 
used for the linear response analysis. 

The classical action for one scalar field (of arbitrary mass and curvature coupling) is 

5^[*,S] = -1 J d^a; V=ff[(V„$)s"''(V6$)+m'$'+ei?$'] , 

where Va denotes the covariant derivative for the metric gab, ^ is the dimensionless curvature 

coupling, and R is the scalar curvature. The path integral over the free scalar field $ is Gaussian 
and may be computed formally by inspection, with the result, 

J [V^] exp (^5^[*,S]) = exp (^-^Tr In G-\g]^ s exp (^Si^lg]) , 

where 

G-'[g] = ^{-n+m^ + ^R) , 

is the inverse propagator of the scalar field in the background metric gab, and the (generally non- 
local) functional 

Sii'ig] = f TrlnG-^b] ^ j d\' ^ {x'\lnG-\g]\x') , 

may be regarded as the effective action due to the quantum effects of the scalar field in this metric. It 
contains an explicit factor of h and records the quantum effects of the free scalar field in the arbitrary 
curved background spacetime gat- No assumption about the smallness of the metric deviations from 
flat spacetime or any other preferred spacetime has been made. 

The expectation value of the energy-momentum tensor of the quantum matter field in this back- 
ground can be formally obtained by the variation, 



(T„,(.)> = --A=_i_5W[^] 



G[g] = / d X V-fl' / d X jht?-. G[g]{x ,x ) 



V-9 V-g J J 5g°-''{x 

= -ifiV-g T^ab G[g]{x, x) 
where Vab is the second order operator (involving covariant derivatives at x and a;'), 



1>abG[g\{x,x') 



1 1 „ 1 2 



-VaVt - -riabU - -m^ilab + iijlabU - VaVb) G\g\{x,x) 

+ i-S:dt + Ivatrr'XV^V^) G[g]{x,x') 



By Noether's theorem, this (unrenormalized) expectation value of Tab is covariantly conserved, 
provided that the effective action 5*^^' [9] is invariant under general coordinate transformations. 
However, (Tat) is divergent because of the singular nature of the limit x' x in (2.5), which requires 
a careful UV regularization and subtraction procedure consistent with coordinate invariance, before 
a finite renormalized value for its expectation value can be defined [3]. 

In physical terms the UV regularization and renormalization procedure mean that the theory is 
not strictly defined at arbitrarily short time and distance scales. The lack of information about the 
physics at those arbitrarily small scales may be absorbed into a finite number of parameters in the 
effective low energy theory at larger scales. Since the effective Lagrangian and energy-momentum 
tensor have canonical scale dimension n (in n spacetime dimensions), the number of parameters is 
given by the number of local coordinate invariant scalars up to dimension n. In n = 4 dimensions, 
these are the parameters of the Einstein-Hilbert action plus the coefficients of the two independent 
fourth order invariants and CabcdC"''"^'^ , where R is the scalar curvature and Cabcd is the Weyl 
tensor, respectively. Thus, in order to renormalize the theory we require the total low energy effective 
gravitational action, 

S.slg] = Si^ig] + J d\^{R-2K) + \j {aCabcdC^''" + /3R') , 

with arbitrary dimensionless constants a and (3. Renormalization means that G , A, a, and /3 are 
at first bare parameters, which may be chosen to depend on the UV cutoff (introduced to regulate 
the divergences in the one- loop term S^^ [g]) in such a way as to cancel those divergences and render 
the total action, Sefilg], independent of the cutoff. Hence the four parameters of the local geometric 
terms (up to fourth order derivatives of the metric which are a priori independent of ?i) must be 
considered as parameters of the same order as the corresponding divergent terms in S^g [g] , which 
from (2.4) is first order in h. Formally, this may be justified by considering N identical copies of 
the matter field, so that s'^' [g] is replaced by A^S^^ [g] and G~^ , A/G„ , a, and f3 are rescaled by a 
factor of A''. In this way all the terms in eq. (2.7) are now of the same order in A'' as A'" ^ 00. 

This formal rescaling by N is carried out at the level of the generating functional of connected 
p-point vertices, Ses[g], (which are the inverse of p-point Green's functions) rather than the Green's 
functions themselves. Therefore, it has the net effect of resuming the quantum effects contained in 
the one- loop diagrams of the matter field(s) to all orders in the metric gab- The large A^ expansion and 
its relationship to the standard loop expansion have been extensively studied in both theory and 
electrodynamics (both scalar and spinor QED) in flat space [51]. The QED case is most analogous 
to the present discussion with the classical vector potential ^4^ replaced by the metric gab- The 
large N approximation (2.7) is also invariant under changes in the ultraviolet renormalization scale 
(by definition of the UV cutoff dependence of the local counterterms which cancel against those 
of the matter action), and is equivalent to the UV renormalization group (RG) improved one- loop 
approximation. 

It is the large N, RG improved one-loop approximation that is necessary to derive the renormal- 
ized semi-classical equations (1.1) with back-reaction, for only in such a resummed loop expansion 
can the one-loop quantum effects of {Tat} influence the nominally classical background metric gab- 
As mentioned in the previous Section, in the ordinary (unimproved) loop expansion the quantum 
fluctuations of the matter can make at most small corrections to the background metric. The large 
N approximation also preserves the covariance properties of the theory, since it can be derived from 
an invariant action functional (2.7). The divergences in {Tab} are in one-to-one correspondence with 
the local counterterms in the action Seff [ff], whose variations with respect to gab produce, in addition 
to the terms in the classical Einstein equations, the fourth order tensors. 



«i?ai,= ^:r^ / d^x^R" = 2gabDR-2VaVbR + 2RRab-^R\ (2.8a) 
J 2 



Hab= ^ d^X CabcdC'"""' = 'iVVCacbd + 2R'"'Cacbd ■ (2.8b) 



/=5 Sg''^ J 



Hence the variation of the effective action (2.7) gives the equations of motion for the spacetime 
metric (for zero expectation value of the free scalar field 3>), 

a ^^^Hab + P ^^^Hab + {Gab + f^9ab) = {Tab} n , (2-9) 

where {Tab) a is the renormalized expectation value of the energy- momentum tensor of the scalar 
field, and aJl the parameters axe now understood to take finite renormalized values. In order for 
the renormalized parameters to be defined unambiguously, we require that any terms of precisely 
the form of the local geometric tensors on the left hand side of (2.9), specified at an arbitrary but 
fixed renormalization scale n, are removed from the expectation value on the right side of (2.9) by 
an explicit subtraction procedure at that scale /x. A concrete example of this subtraction procedure 
in fiat spacetime is given in Section IV. 

It is worth emphasizing that the UV renormalization of the energy-momentum tensor and the 
covariant form of the equations of motion (2.9) are justified by formal appeal to an underlying 
covariant action principle (2.7), whose variation they are. Although particular regulaxization and 
renormalization procedures, such as non-covariant point splitting or adiabatic subtraction, may 
break explicit covariance, the result must be of the form (2.9), with a covariantly conserved {Tat) „ , 
or the procedure docs not correspond to the addition of local counterterms up to dimension n = 4 
in the effective action, as required by the general principles of renormalization theory. Thus, the 
renormalization of the effective action (2.7) suffices in principle to renormalize the equations of 
motion (2.9) and all of its higher variations, a fact we make use of in the next Section. 

The largo N approximation is equivalent to a Gaussian path integration for the quantum matter 
fields, in which the spacetime metric and gravitational degrees of freedom have been treated as 
c- numbers, coupled only to the expectation value of the energy- momentum tensor through (1.1). 
Since the energy-momentum tensor expectation value can be expressed as a coincidence limit of 
local derivatives of the one- loop matter Green's function G[g]{x,x) in the background metric gab 
through eq. (2.5), it requires solving the differential equation G~^[5] o G[g] = 1, or more explicitly 

{-n+m^ + ^R)G[g]{x,x') = ^^^^J^^ , (2.10) 



concurrently with the semi-classical back-reaction equation (2.9). It is the exact solution of this 
equation without any perturbative re-expansion of Gig], and the resulting self-consistent solution of 
eq. (2.9) for the metric gab, that constitutes the principal non-perturbative RG improved feature of 
the large N limit. 

The equations of motion (2.9), which arc the original cqs. (1.1) modified by the additional terms 
required by the UV renormalization of {Tab}, are fourth order in derivatives of the metric. This 
feature, which is not present in QED (but is a general feature of effective field theories which 
are characterized by derivative expansions), has been the source of much discussion in the litera- 
ture [32,33,46,49,47,48]. As is well known from the general theory of differential equations, if the 
order of the equations is changed by adding higher derivative terms, the solutions of the modified 
equations fall into two classes, viz., those that approach the solutions of the lower order equations 
as the new parameters a, /3 — > 0, and those which become singular in that limit. The latter class of 
solutions arc not present in the lower order theory and correspond physically to solutions which vary 
on Planck length and time scales (in order for the higher derivative terms to be of the same order as 
the lower derivative Einstein terms). There is clearly no experimental basis for taking these solutions 
seriously (since they would predict that even empty fiat space is unstable to arbitrarily short length 
and time scale perturbations) [32,33,46,49]. Instead, the modern framework of effective field theories 
suggests that we should regard the Planck scale as the physical UV cutoff which defines the extreme 
limit of possible validity of semi-classical gravity, and that we should confine our attention to only 
those predictions of the theory which involve length scales £ much greater than the Planck length ipi 
[50]. In this regime, the effects of the higher order local terms in eq. (2.9) are suppressed by at least 
two powers of ipi/£, provided the solutions remain regular in the limit of vanishing renormalized 
coefficients a and (5 of the higher order terms. Wo are interested in this paper primarily in defining 
a validity and stability criterion of the semi-classical approximation at length scales £ S> ^pi, and 
only comment briefly on the Planck scale solutions again in the Discussion. 
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III. CTP, LINEAR RESPONSE AND THE STABILITY CRITERION 



In this Section we present our criterion for the vahdity of the semi-classical approximation which 
relies on a linear response analysis. This analysis makes use of {in\in} expectation values which can 
be realized using the CTP formalism [43]. We begin by reviewing a few details of this formalism 
that are needed to derive the causal linear response equation. 

The desired (in\in) expectation values are obtained by integrating the path integral (2.2) along a 
contour from the initial time up to a late time in the future, and then backwards to the initial time. 
This results in a doubling of the field variables with a new CTP index (denoted by capital roman 
letters A,B, . . .), running over the values A = 1,2, that specify the forward or backward part of the 
contour, respectively. After all manipulations are performed the resulting expressions are evaluated 
by equating field variables on the two contours. 

In the CTP formalism the leading order effective action for the matter fields is formally identical 
to eq. (2.4) with the replacement, 

G-'l9]^{g-'lg])^^^CABG-'[g^], 

with 0~^[g'] a 2 X 2 matrix in the CTP indices and cab = diag(l, —1) the CTP metric. Thus, Sf/b] 

depends only on fields of type 1, while G^2ls] depends only on fields of type 2. The signs in the 
CTP metric, cab, keep track of the direction of the time contour, positively directed forward in time 
for field variables of the first type, and negatively directed backward in time for field variables of 
the second type. The corresponding CTP effective action will be denoted by S^^^ . Performing the 
variation of S^^^ with respect to the first CTP component of the metric variable, g^i,, is formally 
identical to that of S^g^fg] in (2.5), and gives {in\Tab{o:)\in) , the unrcnormalized diagonal matrix 
element of the energy- momentum tensor, which is real for hermitian Tat- 

The linear response equation can be obtained by expanding the CTP effective action, in a func- 
tional Taylor series to one higher order, around a given semi-classical geometry gab that solves 
eq. (2.9). Writing 

Sab gab + hab , 

one finds that to second order in hab the CTP effective action is 



oCTPr , ,1 oCTPr 1 , /" j4 SS^ff^lg] , A / N 

SeH lg + h]= S^s [g]+ d X \' Kb{x) 



2 i J 5gf{x)SgB^{x') 



where the first variation vanishes by (2.9). Varying with respect to h^b and then setting — h? — h 
and g^ = g^ = g, gives the linear response equation which is equivalent to the first variation of the 
semi-classical Einstein equations, namely 

^"^^Hab + P ^'^Hab + (Gab + Agab) = 5{Tab) 

= ^Mab^'^hcdix) -\j d^x'^/^^n^f ^=<^(X,X')/Icd(x') , 



where 



is the non-local connected, retarded polarization tensor and Mab'^'' is the purely local part of the 
variation of {in\Tab{x)\in) at x. We follow here the notation of Ref. [6], except for an opposite sign 
convention in the definition of the energy-momentum tensor in eq. (2.2) of that work. 

To demonstrate that (3.5) is indeed the retarded polarization tensor, we carry out the variation 
of the (unrenormalized) CTP effective action for the scalar matter field explicitly, so that 

d{tn Tab{x) in) = = / d X h {x ) 

J Sgf(x)Sg^^{x') 



= ^Mab'"^hcd{x) +iJi J d'^x' {cii'DabSiA[g]{x,x')cAB'D'"'gBi[g]{x' ,x)} hcd{x') 

= jMab'"^hcd{x) +ih J d'^x'y/^ {T>abgii[g]ix,x')'D'"'gii[g]{x',x) - T>abgi2[g]{x,x')T>'"^g2i[g]ix' ,x)} hcd{x') . (3.6) 



The definitions of the various components of the CTP matrix Green's function of the scalar field 

arc [43,51,52] 

Qi2[g]{x,x') = i {in\${x')^{x)\in) = G<{x,x') , 
Q2ilg]{x,x') =i {in\<b{x)^{x')\in) = Gy{x,x) =G<{x',x) , 

gii[g\{x,x) = i {in\T[^{x)^{x)]\in) = e{t,t')G>{x,x') + e{t' ,t)G<{x' ,x) , 

g22[g]ix,x')=g:,[g]{x,x') = ~e{t' ,t)Gy{x' ,x) - e{t,t')G<{x,x') . (3.7) 

Hence suppressing momentarily the spacetime indices in the last line of (3.6), the CTP structure of 
the expression appearing there is 



e{t,t')[G>{x,x')f + e{t' ,t)[G<{x' ,x)f - G<{x,x') G>{x',x) = e{t,t') {[G>{x,x')f - [G<{x',x)f) 

= -\e{t,t'){in\[<^'\x),ii>\x')]\in) , (3.8) 



where we have used the definitions (3.7) and the properties of the Heaviside step function d{t, t') 
for unequal arguments, ignoring possible ambiguities at the coincident points x = x' . Restoring the 
spacetime indices, we find that the non-local term of the variation of {Tab) in (3.4) can be written 

formally as 

-\ j d*x' ^f^'n^:^''"'\x,x')Kd{x) = -y / d*cc' ^e{t,t') {in\[Tab{x),T^''{x')]\in) Kd{x') , (3.9) 

which is real and causal. 

This derivation is still formal because of the singular behavior of the retarded polarization operator 
and possible ambiguities at the coincident points, x = x' . This singular behavior is related to the 
short distance behavior of the formal expressions and their renormalization. The singular behavior 
of commutators of physical currents and their various time ordered products has been recognized for 
some time [44] , and has been discussed in the gravitational context in Ref. [6] . The proper covariant 
definition of the singular functions requires combining the retarded commutator with the first local 
(contact) term, \Mab'^'^hcd{x) in (3.6), in such a way that the divergences in the sum of the two 
quantities can be renormalized via the usual counterterms, namely exactly the same counterterms 
at the level of the effective action which are necessary to renormalize the semi-classical eqs. (2.9) 
themselves. Alternatively, one may calculate the time asymmetric part of the response function, 
which is free of singularities in the limit a; — > k', and define the renormalized time symmetric part 
of the full response function (including the local contact terms) by a covariant regularization and 
renormalization procedure, which gives unique answers up to finite redefinitions of the coefficients 
a and /3 in the fourth order renormalized effective action. It is this latter procedure which we carry 
out explicitly by means of a dispersion integral, after Fourier transforming (3.9) in the flat space 
example provided in the next Section. 

The linearized fluctuation hab(x) obeys an integro-differential equation (3.4) in which the inte- 
gral depends only on the past of x, due to the causal boundary conditions, and which involves the 
two-point correlation function of the matter energy-momentum tensor. According to the general 
principles of linear response analysis, this retarded correlation function is evaluated in the back- 
ground geometry of the leading order solution of the semi-classical equations (2.9). 

The polarization operator, Ill^^'^^'^''{x,x'), is determined by the second variation of the same 
effective action that determines the energy-momentum tensor, and it also obeys the same covariant 
conservation law, 

^-■n'^^-^)-\a;^x') = V,u'^^^'^''''{x,x')=0, for xy^x'. (3.10) 

The equations (3.4) are covariant in form and therefore are non-unique up to linearized coordinate 
(gauge) transformations 

Sgab 5gab + ^aXb + VfcXa , (3.11) 

for any vector fleld Xa- Singular gauge transformations in the initial data for 5gab are excluded, 
and some care is required to decide whether time dependent linearized gauge transformations which 
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grow in time without bound arc allowed or not. Since the action principle is fundamental to the 
present approach, any transformation of the form (3.11), for which the action (2.7) is not invariant 
(due to boundary or surface terms), is not a true invariance and should be excluded from the set of 
allowable gauge transformations of the linear response equations (3.4). 

We now state our stability criterion for the semi-cleissical approximation. A necessary condition 
for the validity of the large A'^ semi-classical equations of motion (2.9) is that the linear response 
equations (3.4) should have no solutions with finite non-singular initial data for which any linearized 
gauge invariant scalar quantity grows without bound. Such a quantity must be constructed only 
from the linearized metric perturbation hat and its derivatives, and it must be invariant under 
allowed gauge transformations of the kind described by eq. (3.11). 

The existence of any solutions to the linear response equations with unbounded growth in time, 
that cannot be removed by an allowed linearized gauge transformation (3.11), implies that the 
influence of the growing gravitational fluctuations on the semi-classical background geometry are 
large, and must be taken into account in the evolution of the background itself. That is to say, if 
the gravitational fluctuations around the background grow, even if they were initially small, then 
the loading order semi-classical equations (2.9), which neglect these fluctuations, must eventually 
break down. 



IV. STABILITY OF FLAT SPACETIME 



Flat spacetime is a solution of the semi-classical Einstein equations for vanishing expectation 
value of Tab and cosmological term, with the quantum matter field in its Lorentz invariant vacuum 
ground state. This is the simplest solution of the semi-classical cqs. (2.9) to which we can apply our 
validity criterion, and for which the polarization operator, U.^^^^^ '^'^ , can be evaluated in closed form. 
In addition to illustrating the application of the criterion to a well-defined specific case, the analysis 
of the normal modes which solve the linear response eq. (3.4), will permit us to roach a definitive 
conclusion on the stability of Minkowski spacetime to quantum perturbations on distance and time 
scales far larger than the Planck scale. 

The linear response equation (3.4) around a Minkowski background, rjab, can be decomposed 
into scalar, vector and tensor components according to the decomposition and projection operators 
defined in Appendix A. The variations of the local tensors appearing on the left hand side of eq. (3.4) 
are given by eqs. (A16). Thus, eq. (3.4) around flat space may be written in the form. 



IGttG, 



SttG, 



hcd{x) 



The non-local vacuum polarization tensor [right hand side of eq. (4.1)] can bo decomposed into ex- 
actly the same two scalar and tensor projections, (see both Appendix A and B), '^''■Ji'^'^) (^''*^) _|_ 
^if "'n'*^ and a Kallen-Lehmann spectral representation [5,4,45] given for the Fourier trans- 
form of each of these two gauge invariant scalar functions, 



n 



(i) (ret) 



f 

J — c 



dw 
2^ 



(a;,fe) 



i = T,S . 



This is a form of Cauchy's theorem for 



icrW(fc°,fe)=7rpW(s)sgn(fc°) 



s = {kr-\k\' 



In Fourier space the non-local real convolution in eq. (4.1) becomes a simple multiplication with 

hcd{k). Its real part is even and its imaginary part odd under time reversal, which is taken into 
account by the sgn(fe'^) function in (4.3), and which is a consequence of the causal, retarded boundary 
conditions of the CTP {in\in) formalism. 

Since the purely local term, M^b"^ , is time reversal invariant, it does not contribute to the imag- 
inary part of the dispersion relation (4.2) of n''"'''''*-'(fc'', fc) which is flnite and well-defined. The 
proper definition of the local term is connected with the renormalization procedure needed to re- 
construct the real part of the Fourier transform of (4.1) from its imaginary part, and in fact, the 
dispersion integral in eq. (4.2) does not exist due to the large k behavior, viz. k^, of the Lorentz 
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invariant spectral functions p'''(s = —k^). This divergent behavior of the unrenormalized disper- 
sion integral (4.2) as oo, is nothing but the ambiguities of the coincident limit x ^ a;' in the 
retarded polarization function in a different guise. The divergent terms are proportional to (5^(a;,a;') 
and up to four derivatives thereof, which by Lorentz invariance must be of exactly the same form 
as the local terms on the left hand side of the linear response eq. (3.4). Thus, these divergences, 
as well as the local term A4^f,'^'^, can be handled by the same renormalization procedure needed to 
define the expectation value of the dimension four operator Tab in (1-1), namely by subtraction of 
the allowed covariant counterterms up to dimension four. In the flat space dispersion integral (4.2), 
this is easily accomplished by subtracting the first three terms in its Taylor series expansion around 
= 0, and defining the renormalized real part of the retarded correlation function by 



Re 



n(')(-')(fc2) 




= Rej 




R- 





t(«) {■■et)/.2N 



= -{ky V 



(ret) , 



d£p«(s) 



QY^ii) (rot) 



dk^ 



^2j-[(i) (ret) 



9(fc2)2 



where V denotes the principal part prescription for the integral when k is timelike (— fe^ = s). The 
subtractions do not affect the time odd imaginary part of the retarded polarization function. The 
integral over s in the real part is now well-defined and UV finite, and may oven be computed in 
terms of elementary functions, in the case of a scalar field of arbitrary mass m > and curvature 
coupling ^. The details of this calculation are given in Appendix B. The three subtractions in (4.4) 
correspond physically to renormalizing the coefficients of the cosmological constant (A), Newton's 
constant (1/G„), and the coefficients of the fourth order terms (a and 0). These are of order (fe^)", 
and {k^Y respectively. The renormalized values of these parameters at fc^ = arc what appear 
then on the left hand side of (2.9). The singular local term, M^i,'^'^ in (4.1) is effectively removed by 
these subtractions in flat space as well, so that the entire linear response eq. becomes well-defined 
and covariant. Indeed had we performed the renormalization at the level of the effective action 
directly then it would be clear that no local term ambiguities appear in the renormalized equations. 

Since the projections onto scalar and tensor modes arc linearly independent (in fact, orthogonal), 
the coefficients of the two projection operators must satisfy the linear response relation separately. 
Transferring the polarization part to the left hand side of (4.1) and taking account of the renormal- 
ization just described, yields two independent dispersion formulae, namely, 



ds 



g {s + k'^ — ie sgn(fc")) 
ds p'^^s) 



{s + — ie sgn(A;'')) 



= 0, 
= 0. 



The two spectral functions are calculated explicitly for the free scalar field with arbitrary mass and 
curvature coupling in Appendix B. For this case, the spectral functions have support only when 
s > 4m^, which corresponds to the two particle threshold for timelike gravitational fluctuations. 
However, some conclusions can be drawn from the two dispersion relations above using only the 
fact that both spectral functions are positive for both transverse tensor and scalar gravitational 
perturbations of flat space. 

Let us examine flrst the tensor dispersion relation. It clearly is always satisfled by k'^ = 0. 
This solution corresponds to the physical, transverse linearized gravitational waves propagating in 
a flat space background. The coefficient of k^ at k^ — is unchanged from the classical value by the 
quantum parameter a and vacuum polarization corrections. Therefore, these linearized gravitational 
waves carry the same energy density in the semi-classical approximation as they do in the classical 
Einstein theory. 

Next we may examine the interior of the brackets to determine if there are any other solutions to 
the tensor linear response equations. Solutions with k^ = — (fe°)^ + > correspond to unstable 
modes with imaginary frequencies, since we can always consider these modes in a frame where fe = 0. 
When k^ > the —ie prescription is not needed, and e may be set to zero. Thus, by making use 
only of the positivity of p^'^\ and assuming that a is non-negative we observe that the bracket is 
strictly positive for k^ > 0, provided 

167raG fc^ + 1 > . 

If a > this is always satisfied, and indeed this constraint on a is required by positivity of the 
energy density —2^'-'^Hoo corresponding to the fourth order CabcdC"''"^''' term in the action. This 



(4.4) 



(4.5a) 
(4.5b) 



(4.6) 
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demonstrates that there are no unstable transverse tensor perturbations of flat spacetime for a > 0. 
Since this conclusion relies only on the positivity of the spectral function p^^\s), it requires only 
causality, a bounded fiauiiltonian, and a well-defined positive Hilbert space norm for the quantum 
matter theory. Hence it is valid much more generally than for the specific scalar field example. 

Going further, we may inquire as to the existence of additional stable solutions characterized by 
propagating tensor wave modes with timelike k {k^ < 0). The bracket in (4.5a) vanishes if 

167raG„fc^ + 1 + S-KG^k^F'-'^^ = , 

where 



Jo 



(s + fe2 — ie sgn(A;°)) 



is a dimensionless function of k^ /m^ (and the sign of fc"), given explicitly for the case of a scalar 
field by (B44). It is clear that if this function remains bounded for all jni^ , the equality (4.7) can 
never be satisfied for 167rG„|A;^| ^ 1, since both the third (polarization) term and first (a) term 
can never be of order unity. In fact, from the explicit form of F^^' for a scalar field, given by eq. 
(B44), 



Re 



9607r2 \ m2 



In I ) as |A;^| — > oo , 



so that the function does grow without bound, but only logarithmically. Hence the relation (4.7) 
cannot be satisfied except at }? approaching , provided > 0. If a < then the preceding 
analyses for fe^ > and fe^ < interchange roles, with the conclusion unchanged. Thus, there are no 
tensor mode solutions of the linear response eq. (3.4) on length scales much larger than the Planck 
length for a massive field theory around flat space, other than the usual linearized gravitational 
waves of the classical theory. On physical grounds one must expect this result to hold for any 
quantum matter field polarization tensor of finite mass obeying the same general properties of our 
scalar field example. 

The logarithmic divergence in the response function f'"^^ when fc^ ^ oo is a consequence of the 
large s (UV) behavior of the spectral function proportional to s^, and is generic, with only the value 
of the finite coefficient of the logarithm dependent on the matter content. However the appearance 
of in the lower limit of the logarithm is a result of our definition of renormalized a parameter at 
k^ = 0, which allows no other scale to appear in the logarithm. This definition is no longer tenable 
in the zero mass limit. Instead, one may renormalize the parameter at an arbitrary non-zero value 
of k"^ = jjL^ , related to our previous definition by 

Then, jx^ replaces in the logarithm of eq. (4.9), and the — > limit may be taken safely by 
maintaining a(/tt^) finite. The conclusions about the absence of solutions to the transverse linear 
response equation (4.5a) at length scales much greater than the Planck length remains unchanged. 

Turning now to the scalar component of the linear response equation, (4.5b), we note the opposite 
sign in the coefficient of the 1 / G^, term, which is the well-known negative metric sign of the con- 
formal factor in the Einstein-Hilbert action. From several different analyses [53,54] it is known that 
there are no physical wavelike scalar excitations of flat space in either the classical or semi-classical 
theory. Physically the reason is that the conformal factor is constrained by the diffeomorphism 
invariance. This implies that the expansion of the gravitational action to second order in the metric 
perturbations about flat space, ft''^' '^^^'^h\^b should be treated as proportional to where x is the 
Fourier transform of a new scalar field variable with no kinetic term in the Einstein-Hilbert action. 
This redefinition may be understood as required also in the covariant path integral treatment of the 
linearized gravitational fluctuations around flat space [54]. The net effect of either the covariant or 
canonical analysis is to remove the overall fc^ = propagating solution from the scalar sector, as it 
is not dynamically allowed by the constraints. The possibility of a finite number of non-propagating 
(global) modes in the scalar sector at A; = cannot be excluded by this argument, and must be 
investigated by imposing appropriate boundary conditions at spatial infinity. We do not treat this 
possibility in this paper. 

Finally, the analysis of the expression within the brackets of (4.5b) shows that there are no scalar 
mode solutions (stable or otherwise) with 4mGj^\k'^\ <C 1, for exactly the same reason as in the 
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tensor case, notwithstanding the sign change in the Einstein term. In the scalar case the expUcit 
form of the response function (B45) yields 



Re 



m 



967r2' 



(1-60'log^ as |A;^|^oo. (4.11) 



which shows that the linear response eq. cannot be satisfied for either sign of unless G„ | \ 
becomes of order unity. This result was first obtained in Ref. [54] for the case of a massive, but 
conformally coupled field with ^ = 1/6. In the conformal case, the large s behavior of the scalar 
spectral function p^^^ is much less severe, as is clear from its explicit form given in eq. (B38), and 
only two subtractions suffice. This corresponds to only a finite renormalization of the jSR^ term in 
the effective action. Ignoring the (3 term completely, the only twice subtracted dispersion formula 
gives then 

in the scalar sector. This form of the linear response equation was used in Ref. [54] to demonstrate 
the stability of flat space to scalar (conformal) fluctuations in the infrared limit, i.e., on wavelengths 
far larger than the Planck length. Indeed, substituting the explicit form of p'*', eq. (B37b) with 
5 = 1/6, into (4.12), shows that it is identical with eq. (4.13) of [54] with 2m*p{s)/3 of that reference 
equal to p^^^ (s) here. 

When ^ 1/6 this argument cannot be used since the scalar spectral function behaves as for 
large s and the twice subtracted dispersion integral in (4.12) diverges. However, analysis of the fully 
subtracted response function behavior in (4.11) above leads to the same result. As in the tensor 
case the logarithmic growth with ]fe^] is generic, with only the coefficient of the logarithm depending 
on the matter content, possibly vanishing in some special cases such as ^ = 1/6. Also as in the 
tensor case the lower limit of the logarithm can be made finite in the limit by redefining 

the renormalized /3 coefficient at a finite fi^ ^ analogously to eq. (4.10). One obtains then the 
same conclusion as for the tensor case, namely that there are no new solutions of the linear response 
equations, stable or unstable, far from the Planck regime, despite the opposite sign of the classical 
Einstein term in the scalar sector. 



V. DISCUSSION 



We have presented a criterion for the validity of the semi-classical approximation for gravity that 
involves solving the linear response equation (3.4), to determine the stability of solutions to the semi- 
classical equations (2.9). If, for a given state and background geometry that solves the semi-classical 
equations, one or more solutions to the linear response equations experience unbounded growth in a 
gauge invariant sense, then the semi-classical approximation is not valid for that particular geometry, 
at least not for that particular state. Clearly this is a necessary, though perhaps not a sufficient 
condition for the validity of the approximation. 

As discussed in the introduction, various methods have been suggested previously to test the 
validity of the semiclassical approximation by making use of the two-point correlation function for the 
energy-momentum tensor. The linear response criterion provides a natural and wcU-dcfiucd way for 
this two-point correlation function to enter into the determination of the validity of the semi-classical 
approximation. Further, linear response involves quantities that lie entirely within the semi-classical 
approximation itself, since the polarization tensor is computed on the semi-classical background 
geometry. The large method, augmented by the causal CTP formulation of the effective action, 
provides a well-defined framework for applying the validity criterion, which is equivalent to a stability 
criterion for the semi-classical solution. 

In the covariant effective action formulation, it is clear that the UV renormalization counterterms 
are the same as those needed to define the semi-classical approximation itself and that there are 
no state-dependent divergences. Although the matter energy-momentum tensor correlator by itself 
suffers from possible ambiguities at coincident points, these are removed by a proper covariant 
regularization and renormalization procedure, which ties these divergences in the correlator at a; = a;' 
to counterterms in the purely local gravitational effective action. The resulting combination of all 
terms in eq. (3.4) then becomes well-defined. 

We have illustrated the use of the stability criterion with the simple example of a quantized scalar 
field with arbitrary mass and curvature coupling in the vacuum state of Minkowski spacetime. In 
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this case it is possible to carry out the analysis to completion and show that flat space is stable in 
the infrared limit. There are no propagating solutions to the linear response eqs. (3.4) around flat 
space, except the usual transverse, traceless gravitational wave excitations of the classical Einstein 
theory, provided we restrict ourselves to solutions with 47rG;y|fc^| ^ 0(1). The exact finite number 
of order unity on the right hand side of this inequality determines the values of /c^ for which new 
Planckian solutions and instabilities will appear. Its value depends on both the matter theory and 
the values of the renormalized coefficients a and /3 of the fourth order terms. The existence of 
such growing modes, which violate the validity criterion proposed here, informs us not that flat 
spacetime is unstable, but only that the quantum fluctuations of the geometry should be included 
in some consistent way at short scales. It is then the semi-classical notion of flat space as a pseudo- 
Riemannian manifold endowed with a smooth metric down to arbitrarily short length scales that is 
breaking down at the Planck scale. The semi-classical approximation which does not incorporate 
the effects of these quantum fluctuations on the mean geometry is certainly not valid in the Planck 
regime. However, the semi-classical approximation contains the signal of its own breakdown at such 
short scales according to the validity criterion, and leads to a completely satisfactory stability of 
flat space for all local perturbations obeying the inequality (4.6). 

That empty flat space with quantum matter in its vacuum ground state should be stable, and 
quantum gravitational effects negligible excluding at the Planck scale, is hardly surprising. It does 
mean that the predictions of the semi-classical approximation at least are not in complete disagree- 
ment with observations in this case. In addition to providing an explicit example of how to handle 
the energy-momentum tensor correlation function by standard renormalization methods to obtain 
well-defined answers, working out this case in detail also provides an important clue as to how the 
validity criterion may fail to be satisfied in more interesting cases. What is required is simply that 
the polarization tensor of the matter fluctuations become singular, i.e., large and unbounded, in 
some region. Only in this way can the natural suppression of G„ fc^ ^ 1 in flat space be overcome. 
A non-trivial example, where new modes may be expected, is a flnite temperature quantum matter 
field in an Einstein-de Sitter model [55]. It would be interesting to apply the validity criterion 
proposed here to a consistent solution of the semi-classical Einstein equations possessing thermal 
matter. 

Further important and interesting examples to which the criterion may be applied are solutions 
with event horizons, such as Schwarzschild and de Sitter spacetimes, as well as more general cosmo- 
logical solutions of the semi-classical equations. If the linearized solutions of (3.4) show any growing 
modes, due to a singular behavior of the polarization tensor and/or the existence of non-trivial 
global modes in such cases, then one would be led to the conclusion that including gravitational 
fluctuations beyond the leading order semi-classical approximation in some form would be required. 
A number of different arguments lead to the conclusion that de Sitter spacetime is not the stable 
ground state of a quantum theory of gravity with a cosmological term [56]. In fact, the two-point cor- 
relation function of the energy- momentum tensor for a scalar fleld was estimated in [6], and argued 
to contribute to a gauge invariant growing mode on the horizon time scale. This proposition could 
be tested by a detailed calculation of the two-point correlation function of the energy-momentum 
tensor and the solutions of the linear response equations (3.4) in de Sitter space. 

A second important application of the criterion is to black hole spacetimes. Ever since the 
discovery of black hole radiance, it has been recognized that the quantum behavior of black holes is 
qualitatively different from the classical analogs at long times, since semi-classical black holes decay 
at late times, while classical black holes axe stable. In the Hartle-Hawking state [57] one can construct 
a static solution to the semi-classical equations (2.9) that is quite close to the classical one near the 
horizon [58-60]. On thermodynamic grounds this state is expected to be unstable [61]. However, 
the stability of this self-consistent solution has not been investigated in a dynamical approach. The 
validity criterion proposed in this paper provides a clear dynamical principle for the stability or 
instability of the self-consistent solutions in both the black hole and de Sitter cases. 
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APPENDIX A: TENSOR DECOMPOSITION AND SPECTRAL REPRESENTATION IN FLAT SPACE 

An arbitrary metric perturbation, hab, axound n dimensional flat space, r]ab, can be decomposed 
in the foUowing way: 

hab = hab + davi + SbV^ + ( dadb - i fjai, +—h, (Al) 

\ n / n 

where is transverse and traceless with respect to the flat metric rjab 

d^hib = = dbhib , and 'n''''hab = , (A2) 

and Va is transverse 

d^vi = . (A3) 

By taking partial derivatives and traces and using these defining properties, the various terms in 
the decomposition can be isolated successively, viz. 

h = v^^hab , (A4) 

w = ^— □ (r;"" - nn-'d'^d'') h^b , (A5) 

vi = U-^{&l-U -^dad') d'^hcd , (A6) 

hab = 



^ TVabTl^'' + ^—a-^dadbd'd'' 



n — 1 n — 1 

1 , 1 . \ n 

hcd , (A7) 



(-d^dbd^ - dldad" + -±-r^,^d^d'' + -±-r^'^''dadb) 
\ n— 1 n— 1 J 



(A9) 



where □ ^ denotes the propagator inverse of □ = rf'^dadb- 
Under an infinitesimal coordinate (gauge) transformation, 

hab hab + daXb + dbXa , (A8) 

the change in hab can be absorbed into a redefinition of the various components of the decomposition 
according to 

h ^h + 2UY , 

w ^ w + 2F , 

Va ^ l^a +Xa , 

, X ■ X 

'^ab ^ i^ab ? 

where Xa has been decomposed into its transverse and longitudinal parts, as 

Xa=X^ + daY , with d^X^ = . (AlO) 

Prom these transformations we observe that the transverse, traceless tensor and the linear 
combination of scalars 

h-nw = ^U''-n-'d^d'')h,a, (All) 

are invariant under infinitesimal coordinate transformations. Hence we may define the projections 
onto the scalar (spin-0) and transverse, traceless tensor (spin-2), gauge invariant terms in the general 
decomposition of the symmetric tensor perturbation hab by 

h^^^ = {r,ab - n-'dadb) (r?=<^ - n-'d^d'') /led 

= i {vab - n-^dadb) ih-nw) = Pif "*/icd , (A12a) 

ftir=C = Pir^'/^cd. (A12b) 

The remaining terms in the decomposition contain all the gauge dependence. We denote the vector 
perturbation [containing both transverse (spin-1) and longitudinal (spin-0) components] by 

ftil' = □ {S'adb + 5tda - U'^dadbd'') d'^hcd = P^V "''hcd ■ (A13) 
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Thus, the general symmetric tensor metric perturbation can be written as the sum of three projected 

components 

i=S,V,T 

The three projectors are orthonormal, i.e., 

p(») e/ p(j) cd _ j-ij p(i) cd 

and complete, and define a unique decomposition (modulo the "^"^^^ conformal Killing vectors in 
flat spacetime). 

Because they are conserved tensors derived from invariant action functionals, all the local tensors 
on the left hand side of the linear response eq. (3.4) must be expressible in terms of only the 
scalar and tensor components of the metric fluctuations. Indeed, by explicit computation in n = 4 
dimensions, 

S'^'^H^b^6n'hif=6n'PZ^^''h.d, 



where hab = Sgab is the metric perturbation (variation). 

The scalar and tensor projectors onto the space of gauge invariant metric perturbations can be 
written in momentum space in the compact forms, 

Pi!^''{k) = ^eab0'"' , 

P^Vik) = I (eiet + eiei) - ^Oabe^' , 



where we have introduced the tensor 9ab, 

kakb 



<'ab = Vab 



which obeys k'^Oab = k 6ab = 0. Therefore, the scalar and tensor projectors are also transverse: 

Because of this property the correlation function of two conserved energy-momentum tensors in 
momentum space 



Kb'^{k) 



^e"'<^-^'\T^b{x)T^\x')) , 



may be expanded in terms of the gauge invariant scalar and tensor projectors only. By conservation 
of Tab this correlator must be transverse, with zero projection onto the vector subspace, i.e., [62] 

Kb^'ik) = Pib'^\k)Il'-'\k) + PZ^ ^'li^^\k) , 

in terms of two scalar functions of momentum k. 

The most convenient way of expressing the retarded correlation function in Fourier space is flrst 
to introduce the spectral function representation [45] in the spin-2 (T) and spin-0 (S) sectors, in 
terms of the Euclidean four-momentum {k^, k) 

— OO 

The retarded correlator is given then by the analytic continuation, ik^ ^ -\- it, i.e., 

n(o (-)(,o^ ^ n«(ifc4 = k° + ie, k) = r ' 

J —OO 
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which corresponds to oq. (4.2) of the text. Moreover, since cr^*^ {lu, k) is an odd function of u (by the 
reahty of H^-*), which otherwise depends only on the Lorentz invariant combination s = oj^ — k ■ k, 
we can define Lorentz invariant positive spectral functions p^'^ (s) by 

(T^'^ (w, fc) = 27V sgn(w) p(*' (s) , (A24) 
and obtain the dispersion formula, 

{l^)-^ab WJ^ s + k^-ie sgn(feO) +^<''' ^"'J^ s + P-iesgn(P) ' ^^^^^ 

which follows by substituting (A24) into (4.2), dividing the integration range over w into positive and 

negative co, renaming the integration variable, and using ds — 2ujduj. Thus, the spectral functions 
p'*'(s) of the two independent scalar functions given in eq. (A21) can be obtained by computing the 
simple correlator in Euclidean momentum and evaluating the imaginary part, after performing the 
specific analytic continuation (4.2), namely 

Im [n^*^ {ik4 = k° + ie)] = 71- sgn(A;°) p''^ {s = -k^) , (A26) 

which is also obtained by continuing the Euclidean 

{k4,f + k-k^k^ -ie sgn(fc°) , (A27) 

to the Lorentzian k^ = rjabk'^k'' = — (fc'')^ + k ■ k. The usefulness of this representation is that the 
imaginary part of the correlator is given by spectral functions which have simple positivity properties 
and which are completely free of ultraviolet divergences. These appear only when the real part of 
the correlator is constructed by the fully covariant dispersion integrals over s in (A25), and may be 
handled by standard methods that make clear their relation to covariant local countertems in the 
effective action. The covariant renormalization of these dispersion integrals by explicit subtractions 
in fiat space is described in Section IV. 



APPENDIX B: GRAVITATIONAL VACUUM POLARIZATION TENSOR IN FLAT SPACE 



The classical energy-momentum tensor for a scalar field in n dimensional Minkowski spacetime is 
given by 

I flat = (1 - 20Va$Vi,f + (2^ - i) 7?„6Vc$V=* - 2e*VaV(,$ + 2^r]ab^VcV'^ - ^VabTu''^^ , (Bl) 



which can be rewritten as 



1 1 1 2 

-dadt - -VabH - -m rjab + S, iVabn -dadb) 



When this is substituted into the Fourier transform of the energy-momentum tensor two-point 
connected correlation function in Euclidean space and the two possible Wick contractions of 
($^(a;)$^(a;'))con are taken into account, we obtain 



na6^'LW= / d"x(T„,(x)T='*(x'))Le 



/ 



ik-{x—x') 



2V^^^{k) V^^'^^ik) H{k) 



{l)cd( 



+ 2Vil\k) V';^lf r'^'ik) + 2VZ^^''' V^'^^^k) h.y{k) 



c' d' 

+ 2Vf^''''' Vff jJ/{k) + 2Vf^^''' Vff KA'^'ik) . 
Here k = {k4, k) the Euclidean momentum, and wo have introduced the following tensors, 

'^abik) = ^{kakb - Vabk'^) + ^{Vabk^ - 2kakb) + ^m^Tjab , 



,T%(2)cd scsd , 1„ cd 

^ab = -^o.h + -^VabV , 



(B2) 



(B3) 

(B4) 
(B5) 



and the notation for the following integrals. 
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H{k) 

r''{k) 



(27r)" p2 ^2 (p _|_ fc)2 _|_ ^2 



y (27r)" p2 _|_ ^2 (p _|_ /j)2 _|_ 



(27r)" p2 + ^2 (p _^ fc)2 + ^2 ' 

.P6 (p= + fc=)(p'' + A;<^) 



We regularize all integrals by means of dimensional regularization in n dimensions [63,64]. By 
introducing a Feynrnan parameter x, the previous integrals are evaluated to yield [64] 



n/2 , 

H{k) = -^—T[2- 

^ ' (27r)" V 



da; \nx' + }i'x(\ — a;)] ^ ^ , 



= I2!1t (2-- 
^ ' (27r)" V 2 



da; + k^x{l — a;)] ^ ^ 



/ {3)cd . (4)cd.,2 . ,2^ 



^a6='W = (2 - I) y dx[m'+k\{l-x)]9-' 



iril''" + rif'')[m' + k'x{l - x)f - -i—r^f "^[m^ + k'x(l - x)] 



n{n - 2) 

X{l-X)j4)cd,2 . ,2 



(n-2) 



n2 (5)cd 



' (n-2) ["^ + ^ - a;)] + a; (1 - x) t^^- j . 

These expressions are given in terms of the following five basis tensors, 

(l)cd / 7 \ cd 
Tab (k) = rjabV , 

'^ab'^'ik) = riabk^k"^ + ■!]'"' kakt , 



'''If"'' (k) = Sakbk"^ + Sakbk" + S^kak'' + S^kak" 

"^ab (^) ~ kakbk k 



In order to obtain the Euclidean polarization tensor we need to compute the corresponding tensor 
products specified in equation (B3). Once we do this we can write the polarization operator (in n 
dimensions and for ^ = 0) as 



(B6) 

(B7) 
(B8) 
(B9) 



(BIO) 

(Bll) 



(B12) 



(B13) 



(B14) 
(B15) 
(B16) 
(B17) 
(B18) 



n„6=<^|^(fe,e = 0)=^F,(/c)rif''(fc) 



where 

1/2 

(27rj" 



^'^^^^Wr^y^-^J I dx[m' + k'x{l-x)]^-' 



- 2n - 4 , 
2n(n - 2) 

k^x 



rn' + k''x{l — x)Y + ~ "'^) + — ) ["''^ + A:^a;(l — a;)] 

feV -A . I f k^ 



+ ^ ((1 - .)^ + .^) + - + ^ -T + - 
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[m^ + A;^a;(l — a;)]^ , 



(B19) 



(B20) 
(B21) 
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Fi{s)- sFi{s) = , 
F2{s)-sF4{s) = 0, 
F3(s)+2F4(s)-sF5(s) = 0. 

In fact the combinations, 2F2 and 3-Fi + 2F2, yield the tensor and scalar spectral functions of 
eq. (A25), respectively. 

Likewise, after computing the x integrals, the ^ dependent polarization tensor becomes 



^2 21 



which is explicitly transverse and proportional to the scalar projector P^^^'^''. Combining this ^ 
dependent contribution with the previous ^ independent part, and recalling (A17) and (A26), we 
may now identify the two independent tensor and scalar spectral functions. 



2. I . . r -12 



> 



Both spectral functions arc positive, as they must be, and agree with results (for ^ = 0) reported 
in [28], and (for arbitrary rn and 5) reported in [29] 

In the case that the curvature coupling takes its conformal value, ^ = 1/6, the scalar spectral 
function does not have terms proportional to or to m^s, and becomes 



2\ -I 



1 - 



€=1/6' ' ' '247r2 V s ' 

which agrees with [54], after account is taken of a relative factor of 2m*/3 in the definition of the 

spectral function p''^-'(s) here, relative to p(s) of that work. 

Finally the integrals appearing in the Kallcn-Lchmann represent ations (A25) are all of the form. 



In,l 



ds 



4„2 S'+1(S + P) 



for fc^ > and n and I integers. By making the change of variables s = 4m^/(l — u^), all integrals 
of this kind may be reduced to linear combinations of 



I„{z) = I„,i=o = 2 / dw 

where 



=0 = 2 / 
Jo 



„2„+2 



4m2 



For / = the Iniz) functions obey the recursion formula, 

2 

^""^^^ " ~ 2n + l ^^-'^"-i(^) ' 



with 



Io{z) = -2 + z log (^^) for z>l 
= -2 + f(^^^ for e>0. 
Using these relations, the response function for the tensor fluctuations can be written as 
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ds 



9607r2 
1 



z'^ + z*Iq{z) 



2 2 

"5 ~ 3" 

46 56 m? 32m^ 



15 3 A;2 (A;2)2 ^ I ^ ^ fc^ 



Am' 



9607r2 

while the corresponding response function for the scalar fluctuations is 



ds 



4m2 



S3(s + fe2) 



1 1 2,^ 2m2 
96^^^ 15 + 3('-^^)-3fc^ + 



-2 + /( 4 
' m2 



As A;^ — » 0, 2: — » 00, and the function Jo (or /) is analytic at 2; ^ = 0. However, as changes sign, 
becomes pure imaginary and 



21 1^-1] tan-i 



4m^ 



for < s = —k^ < Aw? , 



which remains real in this range. Finally when s = —k^ > Am^ , f develops an imaginary part, 
/ ^~2^ ~ ^ ~ ^"^'^ sgn(fc'') for z = sj^ ~ , s = —k^ > Am^ , < 2; < 1 . 



(B44) 



(B45) 



(B46) 



(B47) 



[1] B.-L. Hu, Physica A 158, 399 (1989). 

[2] J. F. Donoghue, Phys. Rev. D 50, 3874 (1994). 

[3] N. D. Birrcll and P. C. W. Davies, Quantum Fields in Curved Space, Cambridge University Press (Cambridge), 1982, and 
references therein. 

[4] See e.g., A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill (New York), 1971. 
[5] See e.g., J. I. Kapusta, Finite Temperature Field Theory, Cambridge University Press (Cambridge), 1989; 

M. Le Bellac, Thermal Field Theory, Cambridge University Press (Cambridge), 1996. 
[6] E. Mottola, Phys. Rev. D 33, 2136 (1986). 

[7] B. S. DeWitt, in Les Houches 1985, Proceedings, Architecture Fundamental Interactions at Short Distances, Vol. 2, 1023- 
1057. 

[8] B. S. DeWitt and C. Mohna. Pan's, Mod. Phys. Lett. A 13, 2475 (1998). 

[9] See however V. Mukhanov, H. Feldman and R. Brandenberger, Phys. Rep. 215, 203 (1992) for developments in a cosmo- 
logical context. 
[10] E. Tomboulis, Phys. Lett. B 70, 361 (1977). 

[11] Our sign and curvature conventions are those of C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, 

San Francisco, 1973). 
[12] A. Ashtokar, Phys. Rev. Lett. 77, 4864 (1996). 
[13] C. Beetle, Adv. Theor. Math. Phys. 2, 471 (1998). 
[14] L. H. Ford, Ann. Phys. (N.Y.) 144, 238 (1982). 
[15] C.-I. Kuo and L. H. Ford, Phys. Rev. D 47, 4510 (1993). 
[16] C.-H. Wu and L. H. Ford, Phys. Rev. D 60, 104013 (1999). 

[17] C.-H. Wu and L. H. Ford, Phys. Rev. D 64, 045010 (2001); e-print gr-qc/0102063. 

[18] L. H. Ford and C.-H. Wu, e-print gr-qc/0102063. 

[19] N. G. Phillips and B.-L. Hu, Phys. Rev. D 55, 6123 (1997). 

[20] B.-L. Hu and N. G. Phillips, Int. J. Theor. Phys. 39, 1817 (2000); 

N. G. Phillips and B. L. Hu, Phys. Rev. D 62, 084017 (2000). 
[21] E. Calzetta and B.-L. Hu, Phys. Rev. D 49, 6636 (1994). 
[22] B.-L. Hu and A. Matacz, Phys. Rev. D 51, 1577 (1995). 



21 



[23] B.-L. Hu and S. Sinha, Phys. Rev. D 51, 1587 (1995). 

[24] R. Martm and E. Vcrdagucr, Phys. Lett. B 465, 113 (1999). 

[25] R. Martm and E. Vcrdagucr, Phys. Rev. D 61, 084008 (1999). 

[26] A. Campos and B.-L. Hu, Phys. Rev. D 58, 125021 (1998). 

[27] R. D. Carlitz and R. S. Willey, Phys. Rev. D 36, 2327 (1987). 

[28] M.-T. Jackol and S. Roynaud, Annalon Phys. 4, 68 (1995). 

[29] R. Martm and E. Verdaguer, Phys. Rev. D 61, 124024 (2000). 

[30] A. Roura and E. Verdaguer, Int. J. Theor. Phys. 38, 3123 (1999). 

[31] D. M. Capper, M. J. Duff, and L. Halpern, Phys. Rev. D 10, 461 (1974). 

[32] G. T. Horowitz, Phys. Rev. D 21, 1445 (1980). 

[33] R. D. Jordan, Phys. Rev. D 36, 3593 (1987). 

[34] J. B. Hartle and B.-L. Hu, Phys. Rev. D 20, 1772 (1979); tbid. 21, 2756 (1979); J. B. Hartle, Phys. Rev. D 22, 2091 (1980). 

[35] G. T. Horowitz and R. M. Wald, Phys. Rev. D 21, 1462 (1980). 

[36] A. Campos and E. Verdaguer, Phys. Rev. D 53, 1927 (1996). 

[37] A. Campos and E. Verdaguer, Int. J. Theor. Phys. 36, 2525 (1997). 

[38] A. Campos and B.-L. Hu, Int. J. Theor. Phys. 38, 1253 (1999). 

[39] E. Calzetta, A. Campos, and E. Verdaguer, Phys. Rev. D 56, 2163 (1997). 

[40] E. Calzetta and E. Verdaguer, Phys. Rev. D 59, 083513 (1999). 

[41] F. C. Lombardo and F. D. Mazzitelh, Phys. Rev. D 55, 3889 (1997). 

[42] P. C. Martin and J. Schwingcr, Phys. Rev. 115, 1342 (1959). 

[43] J. Schwinger, J. Math. Phys. 2, 407 (1961); 

L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964) [Sov. Phys. JETP 20, 1018 (1965)]; 

K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 
[44] See e.g., C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill (New York), 1980. 
[45] G. Kallon, Quantum Electrodynamics, Springer- Verlag (Berhn), 1972. 
[46] W.-M. Suen, Phys. Rev. Lett. 62, 2217 (1989); Phys. Rev. D 40, 315 (1989). 
[47] J. Z. Simon, Phys. Rev. D 41, 3720 (1990); ibid. 43, 3308 (1991). 
[48] L. Parker and J. Z. Simon, Phys. Rev. D 47, 1339 (1993). 
[49] E. E. Flanagan and R. M. Wald, Phys. Rev. D 54, 6233 (1996). 
[50] J. F. Donoghue, Phys. Rev. D 50, 3874 (1994). 

[51] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz, and P. R. Anderson, Phys. Rev. D 50, 2848 (1994); 

F. Cooper, S. Habib, Y. Kluger, and E. Mottola, Phys. Rev. D 55, 6471 (1997). 
[52] R. D. Jordan, Phys. Rev. D 33, 444 (1986). 

[53] J. B. Hartle and K. Schleich, in Quantum Field Theory and Quantum Statistics, T. A. Batalin, C. J. Isham, and G. A. 

Vilkovisky, eds., Hilger (Bristol), 1988; 

K. Schleich, Phys. Rev. D 36, 2342 (1987) 
[54] P. O. Mazur and E. Mottola, Nucl. Phys. B 341, 187 (1990). 

[55] H. Nachbagauer, A. K. Rebhan, and D. J. Schwarz, Phys. Rev. D 53, 5468 (1996). 
[56] E. Mottola, Phys. Rev. D 31, 754 (1985); ibid. 33, 1616 (1986); 

P. O. Mazur and E. Mottola, Nucl. Phys. B 278, 694 (1986); 

I. Antoniadis and E. Mottola, J. Math. Phys. 32, 1037 (1991); 

E. Mottola, J. Math. Phys. 36, 2470 (1995). 
[57] J. B. Hartle and S. W. Hawking, Phys. Rev. D 13, 2188 (1976). 
[58] J. W. York, Jr., Phys. Rev. D 31, 775 (1985). 

[59] D. Hochberg, T. W. Kephart, and J.W. York, Jr., Phys. Rev. D 48, 479 (1993). 

[60] P. R. Anderson, W. A. Hiscock, J. Whitesell, and J. W. York Jr., Phys. Rev. D 50, 6427 (1994). 

[61] S. Hawking, Phys. Rev. D 13, 191 (1976). 

[62] D. G. Boulware and S. Deser, Jour. Math. Phys. 8, 1468 (1967); 

L. S. Brown and A. Zee, Jour. Math. Phys. 24, 1822 (1983). 
[63] See e.g., J. C. Collins, Renormalization, Cambridge University Press (Cambridge), 1984. 
[64] See e.g., B. De Wit and J. Smith, Field Theory in Particle Physics, North-Holland (Amsterdam), 1986. 



22 



